Increase with energy of parton transverse momenta in the 
fragmentation region in DIS and related phenomena. 



B. Blok 

Department of Physics, Technion — Israel Institute of Technology, 32000 Haifa, Israel 

L. Frankfurt 

On School of Physics and Astronomy, Raymond and Beverly Sackler Faculty of Exact Sciences, 

^ Tel Aviv University, 69978 Tel Aviv, Israel 

(N 

o 



Oh 



M.Strikman 



Physics Department, Penn State University, University Park, PA, USA 

Abstract 

The dipole and the DGLAP approximations are combined with the kt factorization theorem to 
i-C demonstrate the fundamental property of pQCD: smaller is the size of the colorless quark- gluon 

configurations in the fragmentation region, more rapid is the increase of its interaction with the 
target as a function of energy. First, we consider two closely related properties of the wave function 
4^ of the projectile (i) the transverse momenta of the quark (antiquark) within the qq pair, produced 

in the fragmentation region by the strongly virtual photon, increase with the decrease of x for fixed 
Q^, (ii) increase of the relative contribution of pQCD to the structure functions as compared to 
^ soft QCD contribution at central impact parameters due to a rapid increase with energy of the 

lO cross section of interaction of small dipoles. Practical consequences of these effects arc presented 

for the cases of the cross sections of DIS and DVCS. We predict that the ratio of DVCS to DIS 
amplitudes should very slowly approach one from above at very large collision energies. Second, 
O we study a closely related phenomenon of the increase of the transverse momenta with the energy 

^ of the characteristic transverse momenta of the gluon/quark configurations responsible for the 

^ transition to the black disk regime. We discuss the impact of this phenomenona on the slowing 

• ^ of the dependence on the initial energy of the coherence length. We demonstrate that a rapid 

projectile has the biconcave shape, which is different from the expectations of the preQCD parton 
c3 model where a fast hadron has a pancake shape. We show that the increase of the transverse 

momenta leads to a new expression for the total cross section of a DIS scattering at very large 
energies, relevant to LHeC and LHC. We discuss the impact of the discovered phenomena on the 
hard processes in pp collisions, and on the dominance of different phases of chiral and conformal 
symmetries in the central and peripheral pp, pA, and AA collisions. 
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I. INTRODUCTION. 



The leading order dipole approximation for the high energy processes in QCD developed 
in Ref. P, cf. also Refs. [2l 13, HI El E] , is the generalization of the pQCD improved parton 
model to the target rest frame description. The leading order dipole approximation provides 
the solution of the equations of QCD in the kinematics of fixed and small, but not too small 
X = Q'^/v and not very small . 

The main aim of the current paper is the study of the properties of the DIS processes in 
the fragmentation region. In particular, we study two closely related properties of the wave 
function of the projectile. 

First, we consider two properties of the perturbative QCD (pQCD) wave function of the 
projectile: (i) we find the increase with energy (1/x) of the transverse momenta (fct) of the 
quark(antiquark) within the qq pair produced in the fragmentation region by the strongly 
virtual photon, (ii) we see the increase of the relative contribution of pQCD to the structure 
functions as compared to the soft QCD contribution at central impact parameters due to 
a rapid increase with energy of the cross section of interaction of small dipoles. We show 
that at sufficiently high energies the transverse momenta kt become larger than As 
a result a new pQCD regime different from the conventional DGLAP [7j, appears at these 
energies. The reason is that the ordering in transverse momenta that led to large \og{Q'^ / h?) 
is not valid any more. Practical consequences of these effects are presented for the cases of 
the cross sections of DIS and DVCS. We predict that the ratio of DVCS to DIS amplitudes 
should decrease with energy and very slowly approach one from above at very large collision 
energies. In addition the ratio cxi/crt oc Q'^/{Ak'^) will also slowly decrease with energy. 
Technically these effects follow from the more rapid increase with the energy of the pQCD 
interaction for smaller dipole and the kt factorization theorem [8l|9]. 

Second, we study a closely related phenomenon of the increase with the energy of the 
characteristic transverse momenta of the gluon/quark configurations in the wave function 
of the projectile responsible for the transition to the black disk/Froissart (BDR) regime. 
We discuss the impact of this phenomena on the slowing of the rate of increase of the 
coherence length as a function of energy. The most striking implication of our analysis is the 
change of the space structure of the wave packet describing a rapid hadron. In the classical 
multiperipheral picture of Gribov a hadron has a shape of a pancake of the longitudinal 
size 1/ IX (where /x is the scale of soft QCD) which does not depend on the incident energy 
[lOj . On the contrary, we find the biconcave shape for the rapid hadron in pQCD with the 
minimal longitudinal length (that corresponds to small impact parameter h) decreasing with 
increase of energy and being smaller for the nuclei than for the nucleons. 

The known rate of increase of the black disk regime momenta leads to the generalization 
of the old Gribov formula for the total DIS cross section in the black disk regime. We are 
able to establish an overall coefficient in the cross section dependence on the log^(l/x) in the 
black disk regime. We determine the relative contributions of pQCD and black disk regime 
into the total small dipole DIS cross section at high energies. This new expression for the 
total cross section of a DIS scattering at very large energies is relevant to LHeC and LHC. 

It is worth noting that the effect of the increase of the transverse momenta in the frag- 
mentation region in pQCD found in this paper is very different from the seemingly similar 
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effect found in the leading aslog{xo/x) BFKL approximation [TT] : for the central rapidity 
kinematics log^ (kf / k^^) oc log(s/so) • The latter is the property of the radiation within a 
ladder, i.e. of a diffusion in the space of the transverse momenta [TT]. Indeed, it has been 
known for some time already that if we look at characteristic transverse momenta in a rung 
with a fixed number N in BFKL ladder, than in the multiregge kinematics the transverse 
momenta do not depend on energy. This fact follows from the derivation of Lipatov diffusion 
equation, where log{kf/kfQ) oc A^-the number of the rung under study. The Lipatov diffusion 
arises since a number of rungs in the ladder increases with the rapidity Y. An alternative 
proof that the transverse momenta do not rise in multiregge kinematics with a fixed number 
of rungs has been given in Ref. [12]. On the other hand the property we are dealing here 
with is the value of the transverse momenta in the wave function of the projectile. 

The rapid increase with the energy of the characteristic transverse scales in the fragmen- 
tation region in the black disk regime has been found first in Refs. [T3| [TU [T51 HHj. The 
prediction of the increase with energy of the transverse momenta in the impact factor , in 
the kinematical domain where methods of the pQCD are still applicable can be considered 
as a precursor of the black disk regime indicating the possibility of the smooth matching 
between the pQCD and black disk regimes. 

The characteristic feature of the LO dipole solution is the approximate Bjorken scaling 
for the structure functions of DIS, i.e. the two dimensional conformal invariance for the 
moments of the structure functions. In this approximation as well as within the leading 
\og{xo/x) approximation, the transverse momenta of quarks within the dipole produced by 
the local electroweak current are restricted by the virtuality of the external field: 

A' < k^t < QV4. (LI) 

Here A = Aqcd = 300 Mev is a QCD scale. It follows from the QCD factorization theorem 
proved in Refs. [HI [5] that within this kinematical range the smaller transverse size d of the 
configuration (the transverse distance between the constituents of the dipole) corresponds 
to a more rapid increase of its interaction with the collision energy: 

a = as{c/S)F^^SxG{x, cjS), (L2) 

here F"^ = 4/3 or 9/4 depends whether the dipole consists of color triplet or color octet 
constituents and is a SU(3) Casimir operator, G is an integrated gluon distribution function 
and c is a parameter c = 4 -^ 9. It is well known in the DGLAP approximation that the 
structure function G{x,Q'^) increases more rapidly with 1/x at larger Q^. This property 
agrees well with the recent HERA data. 

Our main result is that the median transverse momenta k^ and invariant masses of the 
leading qq pair in the fragmentation region in the pQCD regime grow as 

k^ ~ am/{xo/x)'^^"\ 

~ 6(g')/(xo/x)^M{Q^), 

(L3) 

Here kf and are the median squared transverse momentum and invariant mass of the 
quark- ant iquark pair in the fragmentation region. (The median means that the configu- 
rations with the momentum/masses less than the median one contribute half of the total 
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crosssection). The exponential factors A and Am are both approximately ~ 0.1. These fac- 
tors are weakly dependent on the external virtuality Q^. Their exact values also depend 
on the details of the process, i.e. whether we consider the DIS process with longitudinal 
or transverse photons, as well as on the model and approximation used. This dependence, 
however, turns out to be very mild. The exact form of A(Q^), and Am(Q^) is given below. 
The function hiQ"^) ~ a((5^)/4 for longitudinal photons due to relation 

z{i-zy ^ ' 

and the dominance of the z = 1/2 configurations. For transverse photons however a{Q'^) <^ 
b{Q^) for realistic energies due a dominant contributions of configurations with 2; ~ 0, 1 in 



the total crosssections, i.e. ^ Ak'f due to the relation 1.4 {z is the fraction of the total 



momentum of the dipole carried by one of its constituents). The characteristic invariant 
masses of transverse qq pair are of the order at realistic energies. 
We carry out the calculations in a wide region of a phase space: 

10"^ <x< 10-^ 3 GeV^ < < 100 GeVl (1.5) 

We also obtain the increase of the median transverse momenta and invariant masses for the 
jet distributions for both longitudinal and transverse photons. 

Although the increase of the transverse momenta in the fragmentation region was ob- 
tained in this paper in the LO approximation of pQCD it should be valid in the NLO as 
well. Indeed, it is known that the NLO radiative corrections to the impact factor are small 
[171 [18]. There are indeed large radiative corrections to the ladder however these corrections 
are already included in the LO CTEQ functions, that were matched to phenomenological 
data for the experimentally studied region of Q^,x. 

Let us note that the transverse momenta for the case of the transverse photons were 
calculated in the framework of the perturbative QCD. In order to ensure separation from 
soft physics (the AJM model), we calculated the average parton momenta with the cutoff 

~ 0.35 GeV^. The relative contribution of the AJ quickly decreases with energy, although 
it is still considerable at HERA energies and few GeV^. 

Finally, we discuss the possible applications of our results to pp, pA collisions at the LHC. 
We explain that characteristic phenomena at central impact parameters are the blackening of 
the interactions, the increase with the energy of the parton momenta, and suppression of soft 
QCD contributions. At sufficiently large energies the collisions at central impact parameters 
will be dominated by a phase with unbroken chiral and conformal symmetries . On the 
contrary, the peripheral collisions will be dominated by a phase with the spontaneously 
broken symmetries. Applicability of this picture to QCD physics at LHC will be subject of 
further research. 

The paper is organized in the following way. First, in section 2 we derive the spectral rep- 
resentation for the total crosssection of the DIS. The derivation is done within the framework 
of dipole model and kt factorization and generalizes LO DGLAP [7] and BFKL [H] approx- 



imations. Next in section 3 we derive eq. 1.3 in double log approximation. This derivation 



has an advantage that the increase of the transverse momenta can be seen analytically, in 



a model independent way. Next, in section 4 we derive eq. |1.3| in LO approximation using 
CTEQ5 and CTEQ6 structure functions. We find the kinematical boundary of the new 
pQCD regime. We also discuss the influence of the soft (AJ) contributions on our results. 
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In section 5 we applied our results to DCVS processes and study the dependence of the 
cross-section of these processes on and energy xb- In section 6 we derive the generali- 
sation of Gribov formulae- the formula for total DIS cross-section in the black disk region, 
that takes into account the rapid rise of the transverse momenta in the fragmentation region 
in the black disk regime. In section 7 we study the coherence legth and the space-time form 
of the nucleons/nuclei in the black disk regime. In section 8 we discuss the experimental 
consequences of our results. 

Some properties of the dipole approximation in QCD and DIS we use are discussed in 
appendixes. In appendix A we remind the reader the basic ideas of the derivation of the 
dipole model, in particular we give a detailed derivation of the momentum representation 
of the total crosssection. In appendix B we remind the reader the basic formulae of the 
AJM model. In appendix C we recalculate the characteristic transverse momenta for the 
black disk regime, using CTEQ6 distribution functions and two-gluon nucleon formfactor 
parametrization from ref. [131 IB]- In appendix D we present the parametrizations of the 
structure functions that we use in this paper. In appendix E we present some interpolation 
formulae that describe the data on median momenta as function of x and as described 
in tables 1,2. 



II. THE TARGET REST FRAME DESCRIPTION. 



In this section we shall use the dipole approximation to derive the momentum represen- 
tation for the total cross-sections of the DIS for both transverse and longitudinal photons. 
Within the LO approximation the QCD factorization theorem allows to express the to- 
tal cross section of the scattering of the longitudinally polarized photon with virtuality 

^ ^QCD off ^ hadron target as the convolution of the square of the virtual photon wave 
function calculated in the dipole approximation and the cross section of the dipole scattering 
off a hadron [H |5l [191 120] • In the target rest frame the cross section for the scattering of 
longitudinally polarized photon has the form (see appendix A in this paper for the detailed 
derivation of this formula): 



a(72 + T^X) 



127r2 



d'^ktdz(^lj^*^{kt,z) a{s,kt) ip^*{kt,z)y (2.1) 



Here a is the dipole cross section operator: 



a = F^- 7T^a,{4k^){-At) ■ xG{x = (M^ + Q^)/sAK 



(2.2) 



At is the two dimensional Laplace operator in the space of the transverse momenta, and 
is the invariant mass squared of the dipole. In the coordinate representation a is just a 
number function, and not a differential operator as in the momentum representation. 

In the leading log(a;o/a;) approximation a similar equation arises where the cross section 
is expressed in terms of convolution of impact factor and unintegrated gluon density. In 
practice, both equations should give close results. Integrating by parts over kt it is easy to 



rewrite Eq. 2.1 with the LO accuracy in the form where the integrand is explicitly positive: 



a(7l + T^X) = ^J «s(4fc,2)c?'M^(VV^^.(A;i,^)|/(s,z,fc2)|v^,*(fci,z)), (2.3) 



here 



f = {An'/3)a,{Akt)xG{xAkt) 



(2.4) 
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In the derivation we use the boundary conditions that follow from the fact that the photon 
wave function decreases rapidly in the kf ^ oo limit and that the contribution of small kt 
is the higher twist effect (see Appendix A for details). 

The cross section of the interaction of the longitudinal photon can be rewritten in the 
form of spectral representation by explicitly differentiating the photon wave function : 

a.(x,g^) = 6vr ""-y^ ^ J dM^ J dza^iAM^zil - z)) ^^f^ ^^^^ ■ gix, (2.5) 

Eqs |2.3[ |2.5| follow from the kt factorization theorem for high energy processes. Here 



g{x,Q'^) = xG{x,Q'^).Tlms Eq,2.5 is the generalization of the DGLAP approximation to 
the domain of high energies where the processes with large and increasing with energy 
play an important role. 



Note that in Eq. 2.5 the dependence on z comes explicitly only in the argument of 
the coupling constant. Practically, for the do mina nt configurations z ~ 1/2 ( and hence 
4fc^ ^ M^). Then the dependence on z in Eq. 2.5 is absent and one can consider Eq. 2.5 



(with z=l/2) as a spectral representation of a electroproduction amplitude for longitudinal 
photons. The spectral representation of the electroproduction amplitude over is a general 
property of a quantum field theory at large energies, where the coherence length significantly 
exceeds the radius of the target T [211 122] . The DGLAP approximation ensures additional 
general property: the smaller size of the configuration in the wave function of the projectile 
photon leads to the smaller interaction with the target, but this interaction increases with 
the energy more rapidly. Such spectral representation is not valid for transverse photons due 
to the AJM effects, since the configurations with smaller z dominate (see below), and the 
integrand in the corresponding equation contains the contribution of different configurations 
(contrary to the preQCD spectral expansions). 



Eq. 2.1 can be equivalently rewritten in terms of integration in kf and z as 

a.(x,Q^) = 6. ""-y^ ^ / ^^X(4fen^^(l-^)^ (^.^^,;(,_^)), -^(^,4A:n. (2.6) 

where x is given by {kt/{{z{l — z)) +Q'^)/s. Here we take into account explicitly the (rather 
weak) 2;-dependence of the integrand. 

In the NLO approximation the structure of formulae should be the same except the 
appearance of the additional qqg, ... components in the wave function of photon due to the 
necessity to take into account the QCD evolution of the photon wave function [TU] . 

The similar derivation can be made for the scattering of transverse photon in configura- 
tions of spatially small size. In this case the contribution of small kt region (Aligned Jet 
Model contribution) is comparable to the pQCD one. The main interest in this paper is 
in the region of high energies (HERA and beyond) i.e. sufficiently small x, and small Q^, 
where pQCD contribution dominates because of the rapid increase of the gluon distribution 
with the decrease of x. The pQCD contribution into the total crosssection initiated by the 
transverse photon has the form: 

(Ty = OTT — 



12 



r ^^/^M^«.(4M^.(l - ^) f!^t~zflM4Q^r ■ - 

(2.7) 
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The invariant mass representation for the transverse photon does not give objective infor- 
mation on the parton transverse momentum distributions. The reason is of course the well 
known fact, that there are configurations with large invariant masses, but small transverse 
momenta k^. These configurations are the perturbative analogues of the AJ configurations. 



It is useful to rewrite Eq. 2.7 in terms of integration over and z. The reason is that the 
median invariant mass and transverse momentum for transverse photons are not connected 
like in the longitudinal photon case, and we need two different equations to determine them. 
Cross section initiated by the transverse photon has the form: 



(Tt = 67r- 



12 



(2.8) 

We include a contribution of the aligned jet configurations by imposing a cut off in transverse 



momenta. In particular, in the numerical calculations using Eq. |2.8| we introduced a cutoff 
in the space transverse momenta kf = M'^z{l — z) > u,u ^ 0.35 GeV^. The contribution of 
smaller k'f in the total crosssection was calculated using the AJM model. 

III. THE DOUBLE LOGARITHMIC APPROXIMATION. 

In this section we analyze the new properties of the pQCD regime within the double log 
approximation. The advantage of this approximation is that it will allow us to perform some 
of the calculations analytically. Other calculations will be made numerically but using the 
expressions that are known analytically. 

In the double logarithmic approximation the structure functions are given by [23] 

xGix,Q') = I dj/i2m)ix/xoy-\QVQlr^^\ (3.1) 
where the anomalous dimension is 

<3 - 1) 

To simplify the calculation we assume, the initial condition for the evolution with Q^: 

G{x,Ql) = 5{x-l). (3.2) 
In the saddle point approximation one finds |23j : 

xG{x, Q') = i^g^^^^^^°3/4 \/4«s(gg)(A^c/vr) log(gVQ^) ^og{xo/x). (3.3) 

Structure function of a hadron is given by the convolution of this kernel with the nonper- 
turbative structure function of a hadron in the normalization point Q"^ = Qq- Note that G 
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is increasing with .This is the pQCD contribution where virtuahties of exchanged gluons 
are large. 

In the analysis of energy dependence of the parton momenta it is legitimate to neglect 
the pre-exponential factor, since the absolute value of G as well as the pre-exponential factor 
weakly influence the transverse scale, and its evolution with energy: 

xG{x, = exp ^Aas{Ql)Nj7: \og{Q^/Ql) log(xo/x). (3.4) 

We shall consider here the double logarithmic approximation without the running coupling 
constant. Taking into account the running coupling constant leads to similar results. We 
shall not present them here, since they are more cumbersome, and are essentially reproduced 
when we shall tackle the full LO calculations in the next section. 



A. Energy dependence of the quark transverse momenta for fragmentation pro- 
cesses initiated by longitudinal photon. 

We shall find analytically the scale of the transverse momenta in the limit where s ^ 
^ . For certainty we restrict ourselves to the contribution of light quarks. 

At large the crosssection for the scattering of the longitudinal photon is dominated 
by the contribution of the spatially small dipoles, so it is legitimate to neglect the quark 
masses. In this limit the cross section is proportional to 

aLOcQ^ J dM^n{M\ s, Q^), (3.5) 



where the function n(M^, s, Q"^) is given by Eqs. 2.5, 2.6 
n(M^.,Q^) = «.(MV4)p^^^^ 

X exp(74a,(Q2)(iV^/,r) log(MVQg) log((xos/(M2 + Q2))^ 

(3.6) 

where Xq = Q'^/sq. Here we keep only large terms depending on (we ignore the 
independent normalization factor irrelevant for the calculations below). 

Let us show that the maximum of n(M^, s, Q"^) increases with the energy. At very high 
energies n is proportional to 

n ~ exp(log«,(MV4) + logiM^Q^) - 41og((Q2 + M^)/Q^) 
+ ^j4as{Njn){\og{s/so) - \og{{Q' + M^)/Ql)) log(MVQ2)). 

(3.7) 

In the limit of fixed but very large energies, log(s/so) ^ log(((5^ + M^) /Qq )). Let us 
assume that for the maximum: ^ Q^. The maximum of the expression 3.7 under this 
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assumption can be found analytically. We shall see that this assumpti on i s indeed self- 
consistent. Indeed, under the latter assumption we can differentiate Eq. 3.7 over and 
find for the maximum: 

M2 = M2(s/so)"=(^^/")/^ (3.8) 
Here Ql ~ and sq ~ Q^- We will refer to this extremum value of as . We see 



from Eq. 3.8 that the condition M( » Q is self-consistent at very high energies. 



At the extremum n oc (as(M^/4)/Mf exp((A^c/7r)(as/3) log(s/so)). Therefore 
daL 



dM 



|M2=Mf ^ a.(M//4)(gVM{')(exp((iV,/7r)(a,(g^)/3) log(s/so)) (3.9) 



However, the position of the maximum of the integrand is not sufficient to characterize the 
relevant transverse scales as a large range of is important in the integrand. In particular, 
calculation of second derivative shows that dispersion over = is large. The width of 

the distribution over \og{M^/M^) is ^2/3 \og{Ml/M^). 

Hence we need to determine range which gives most of the integrand support. For 
certainty, we define the range of < Mf which provides a fixed, say, 50% fraction of the 
total perturbative cross section. Let us estimate how this scale increases with the energy in 
the double log approximation. First, let us consider the total cross section. The upper limit 
u of integration over is determined by the kinematic condition through tmm, giving that 
the allowed invariant masses <^ s. We choose upper limit of integration as 

< M\^, = 0.2s, (3.10) 

from the cross section of diffraction although the result of numerical calculations is insensitive 
to the upper bound because essential are significantly smaller. In fact the integral for 



the cross section converges long before the upper limit of integration 3.10 is reached (see the 
discussion below). 

Let us first calculate the median scale semi analytically. Within the double logarithmic 
approximation, and assuming that the conditions log(s/so) ^ log(((5^ + M^)/ {Q^ + Mq), is 
still valid for the relevant M^, the integral for the cross section can be written as: 

rlog(n/M2) 



a{u) = (Q'/Qt) / dlog(MVQ^)a,(MV4)exp(-21n(MVg, 
Jo 

+ V(4«.Ar,/7r) log(MVg2) log(s/so)). 



(3.11) 



Here u is the upper cut-off in the invariant masses. Introducing the new variable t 
log{M'^/Ql), we obtain: 



{QVM^) r^"^ dtas{tMl/A) exp(-2t + y^(4a,Ar,/7r) log(s/so)t), (3.12) 







where k{u) = logfu/so). The integral for the total cross section is given by the equation 

similar to Eq. 3.12, with the upper integration limit being replaced by n{s) = ^log(0.2s/so). 
Note that the essential scale of integration is determined by exponent, and is very weakly 



influenced by the exact value of a upper cut-off. The integral 3.12 is actually the error 
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function [21], which can be easily evaluated numerically. Requiring that it gives one half of 
the cross section we find 

~ gg^s/so)°-''"='^=/". (3.13) 

Evidently, for sufficiently large s our initial assumption log(M^/QQ) ^ log(Q^/(5o) f^^y 
self-consistent. This is because the decrease of n with due to 1/M^ terms in the integrand 



of Eq. 2^ is partially compensated by the rising exponential, giving a relatively slow decrease 
of n to the right of the maximum of the integrand. 

Note that the rate of the increase of with s is much higher than for due to the 
slow decrease of the integrand with M^. The cross section of jet production with at this 
interval also increases with the energy as 



da 



2 ~ 



s/so)°-''"^^=/". (3.14) 



The analytical calculations supply the pattern for the behavior of the transverse momenta 
in double log approximation. In order to understand the dependence of the median scale on 
both the energy and quantitatively in the double logarithmic approximation we made 
numerical calculation of the characteristic transverse momenta using the DGLAP double 
log structure function. We find that the increase rate of the transverse momenta indeed 
does not depend on the external virtuality Q^. Considering the wide interval of energies and 
s = 10^ 10^^ GeV^, and 20 < < 200 GeV^ we obtain the approximate formulae: 

M^{x) ~ 0.7Q2) exp(0.17((4a,iV,/7r) log(xo/a;))°-^^). (3.15) 

Here Xq = 0.01. We present this estimate only for illustrative purposes, since the double 
logarithmic approximation is semirealistic only. Still our results indicate that for external 
virtualities < 100 GeV^ and energies which can be reached at LHeC the onset of a new 
pQCD regime may take place. 

The cross section of the jet production at this scale also increases with the energy as 

'^'^ 'm^=m? - iGix,Q')/Q'){Mfix))) (3.16) 



dM^ 



where M^{x) is given by Eq. 3.15 The median transverse momenta ~ /4, due to the 



dominance of symmetric configurations. 



B. Transverse photon: the characteristic transverse scale in the photon fragmen- 
tation region. 

The main difference between the longitudinal and transverse structure functions in the 
DIS is the presence of the strongly asymmetrical in z configurations due to the presence of 
the {z{l — z))~^ factor in the spectral density. As a result there is a competition between two 
effects. One is a slower decrease of the spectral function with (by the factor M'^/Q^), 
leading to the more rapid increase of the characteristic transverse momenta for the symmetric 
configurations. Another effect is the presence of the asymmetric {z — >■ 0) configurations 
which are characterized by the small transverse momenta kf for a given invariant mass M^. 
For such configurations the rate of increase of the gluon structure function with energy is 
small. 
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Let us first show that the transverse momenta increase rapidly for symmetric configura- 
tions. The spectral representation for the transverse photon for symmetric configurations is 
proportional to 



(M2 + Q2)4 



X 



Aas{N,/7r){log{s/so) - log{{Q^ + M^)/{Ql)) \og{M^/Ql)). 

(3.17) 

In the high energy limit, when 3> Q^, we find for the dependence of the maximum of n 
on energy: 

Ml ~ QI{s/sq)Y^^'''^'''^I\ (3.18) 

This is twice as fast increase as for the case of the longitudinal photon. increases 
with s at high energies and thus the condition ^ is perfectly self-consistent at very 
high energies. In addition we can calculate the total cross section in the same approximation 
semi-analytically getting the error function and obtain the rate of increase [s / Sq))'^'^^^'^'^''^''/'^\ 
which is twice that for the longitudinal case. 

These two results are applicable to the symmetric configurations only. On the other hand, 
at least at achievable energies, the dominant contributions for transverse photon crosssection 
are asymmetric, with z close to or 1. In order to take these configurations into account we 
performed a numerical calculation using the gluon distribution function within the double 
log approximation. The result is that the characteristic median scale M| increases as 

The value of the exponent is 0.12 for the beginning of the studied energy range s ~ 10^^10^^ 
GeV^, and decreases to 0.09 at the upper end (for typical Us = 0.25). Thus the rate of the 
increase with the energy is approximately the same as for longitudinal photons for not very 
high energies. For very high energies the symmetric configurations win over asymmetric 
ones, leading to the increase of the average transverse momentum squared which is twice 
as large as in the longitudinal case. The precise determination of the scale A1q{Q'^) is 
beyond the accuracy of this approximation. Effectively we obtain the dependence M| ~ 

0.7Q2(a;o/a;)0.1(4a.iV,/^)_ 

One can also estimate the rate of the increase of the jet production cross section : 

dar/dM^, - {xG{x, Q^)/Q%1/3){1 + 0.5{xo/xf-^^). (3.20) 

We found a rapid increase of the jets multiplicity with energy. Thus the rate of the increase 
with energy of the transverse momenta of quarks in the current fragmentation region for 
transversely polarized photon is significantly more rapid. Consequently we find that a^/aT ~ 
asQ'^/M'^ being numerically small should slowly decrease with energy at sufficiently high 
energies . 

We conclude that it is possible to show analytically that for very high (asymptotic) ener- 
gies the relevant invariant masses extend well beyond and increase with the energy. The 
direct numerical calculation of the Mf scale shows that the rate of increase is independent 
of external virtuality. 
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Note that due to the significant contribution of the asymmetric configurations the median 
transverse momenta is much smaller than M| /4. The simple numerical calculations using 

= M'^z{l — z), shows that the average transverse momenta (including nonsymmetric 
configurations) rapidly increases like a{Q'^)/{x/0.01)^'^'^, with the exponent once again is 
independent of external virtuality, and a{Q^) ~ 0.6GeV^ + 0.02Q^, i.e. is much smaller 
that especially for small virtualities. 

IV. THE CHARACTERISTIC TRANSVERSE MOMENTA IN HARD FRAG- 
MENTATION PROCESSES IN LO APPROXIMATION. 

In the previous section we discussed the quantitative dependence of the characteristic 
transverse scales of DIS in the double logarithmic approximation in the dipole model at 
unrealistically high energies, where the analytical calculations are possible. Here we carry 
out the calculations for realistic energies and realistic structure functions. The numerical 
results indicate that the effects discussed above are manifest even at the energies of the order 
s ~ 10^ ^ 10^ GeV2. 

We will also consider the extrapolation of our results to energies of the order s ~ lO^GeV^. 
These energies are unattainable at existing facilities. The proposed LHeC collider may reach 
the invariant energies of order 10^ GeV^ [25]. However these results are interesting from the 
theoretical point of view- probing the limits of the pQCD. 

The challenging and unresolved problem is how to use resummation methods at extremely 
small X [261 l2Z| to evaluate dependence on energy of parton distribution in the current frag- 
mentation region. One can substantiate this point by evaluation of the number of radiated 
gluons in the multiRegge kinematics [IS], obtaining 1-2 emissions at HERA energies. How- 
ever at extremely small x where number of gluon radiations would be sufficiently large and 
therefore these models can become applicable. This interesting problem is beyond the scope 
of this paper. 



A. The longitudinal photons. 

In the case of longitudinal photons we have considered the characteristic median/average 
transverse momenta scale, that corresponds to the half of the total crosssection ai- This 



scale is determined from Eq. |2.5| by first integrating over z for given kt, and then analyzing 
the corresponding jet distribution. The results of the numerical analysis using the CTEQ6 
gluon distribution functions are given in the Table 1. In Figur^we present the characteristic 
graphs for the ratio 

R{kf) = (4.1) 



where (T{k^) corresponds to the result of integration of Eq,2.6 



over transverse momenta 

< kf. We see from Fig. [T| that for fixed kt R{kt) slowly increases with the increase of 
the energy. The results based on using CTEQ5 parametrization are qualitatively similar, 
although the increase of median k^ with the energy is more rapid. The energy depen- 
dence of median kf can be described with a very good accuracy by an approximate formula 
(a;/0.01)°-°^+° °25i°g(QVQg). Here = 10 GeV^xo ~ 0.01. The power increases from ~ 0.04 
at ~ 5 GeV^, to 0.09 at ~ 100 GeV^. The interpolation expression which provides a 
good description of the results presented in Table 1 is given in Appendix E. For CTEQ5 the 



12 



exponent in Eq. El increases to 0.1 at = lOOGeV^ instead of 0.09. This is consistent 
with the enhanced rate of the increase of CTEQ5 structure functions as compared to the 
CTEQ6 ones (see below). 

These results allow us to estimate the scales, where one expects the appearance of the 
new QCD regime, i.e. one has to use the kt factorization approach. Indeed, the DGLAP 
approximation is based on the strong ordering in all rungs of the ladder, in particular in 
the first rung (the impact factor in the 4kt factorization language ) we must have AAqqy) ^ 
< Q^. It is clear, this ordering can not hold, once the median 4kf becomes of order Q^. 
Then we obtain the condition (using CTEQ6 distribution functions): 

4a(Q2)/(a;/o.01)0-o™5iog(QVQg) ^ q2 (4 2) 

Here the function a corresponds to the transverse momenta at x = 0.01. 

The analysis of Table 1 and Figjl] shows that for = 5 GeV^ one gets from eq. 4^ 
X ~ 10~^, for = 10 GeV^ one gets x ~ 10~^, which may be reached at LHeC For 
larger we are however beyond the realistic energies: say for ~ 20 GeV^ we need 
X ~ 10~^. The use of CTEQ5 gives qualitatively the same results (for = 30 GeV^ we 
obtain x ~ 10^^. Thus we may hope to observe the onset of the new regime for the kt 
dependence analyzing small x jet distributions at LHeC/LHC. 



The similar analysis can be made using Eq. 2.5 and considering the invariant masses M 



r2 



We obtain that our results for median invariant masses are ~ (4 5)A;|, confirming the 
dominance of the symmetric configurations. 



B. Transverse photons. 



We perform the numerical analysis for the transverse photons using eqs. 2.7[2.8 in the 



same fashion as for the longitudinal photons. The results for the median transverse momenta 
are given in Table 2. In figure^ we depicted the characteristic function R{kf) given by Eq. 



4.1 that gives the characteristic momenta as a function of x for typical value of = 40 GeV^. 
The characteristic energy dependence in Table 2 for median k'f is (^x /0.01)^'^^'^^''^^^^°^^'^ ^^0) 
where Xq = 0.01, Qq = 10 GeV^. A complete expression is given in appendix E. The curves 
in Fig. [2] clearly show that the characteristic momenta increase with the increase of 1/x, as 
the corresponding curves slowly shift to the right. 

We see from the Table 2 that the average transverse momenta for longitudinal photons 
is significantly larger than for transverse photons. On the other hand, the invariant masses 
for transverse photons are always significantly larger than 4kf. This is due to the large 
contribution of the AJM type configurations with 2; ~ 0, 1 . Since = kt/{z{l — z)), a 
more slow increase of than of k^ is consistent with the slow increase of average z towards 
1/2, i.e. the symmetric configurations become dominant, but only at asymptotically large 
energies. 

Once again, we can estimate the boundary of the region where the direct DGLAP ap- 
proach stops being self-consistent. Assuming kf ~ Q'^/4:, we obtain that the boundary for 

= 3, 5, 10 GeV^ is reached at x ~ 10~^, 10~^, 10~^. For higher this boundary lies at 
unrealistically high energies. The use of the CTEQ5 parametrization gives qualitatively the 
same results. 

So far we considered only perturbative QCD contribution, and the median transverse 
momentum in Table 2 is determined relative to the total perturbative cross section, i.e. the 
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one starting from the cut off m = 0.35 GeV^. However, as we mentioned above, there is 
also important AJ contribution in a total cross section. In this paper we take them into 
account using the AJM model [20] (see also Appendix B). This contribution is given in 
Table 3. It is well known that at low x AJM gives a dominant contribution into cross 
section. For x ~ 0.01, ~ 1 GeV^ we see from Table 3 that AJM gives 70% of the 
total cross section. However even at HERA energies the contribution of AJM into the total 
cross section remains significant. Note that the median at small virtualities at HERA 
energies significantly decreases if we calculate it using the cross section that includes both the 
pQCD and soft (AJM) contributions. For example, at ~ 10 GeV^ the median transverse 
momentum squared decreases by almost a factor of two down to ~ 0.65 GeV^. 



V. DEEPLY VIRTUAL COMPTON SCATTERING. 

As the application of the formulae obtained in this paper we shall consider the DVCS 
processes 7 + p — > 7* + p. We shall show that the slow increase in the median transverse 
momenta leads to the slow decrease of the ratio R = ^dvcs/^cs with energy to the limiting 
value equal one. 



The DCVS amplitude is described in pQCD by the same formula 2.7 as the amplitude 
describing total cross section of DIS at given x,Q^ but with the substitution in Eq,2^ of 
the wave function of virtual photon by wave function of a real photon, i.e. = 0. 

As a result in pQCD R has the form : 

^ /o^ dzJdMWsjM'zjl - z)){l/{M' + Q^) ■ g{x, M'). 
""^^^ S^dzSdM^as{M^z{l-z)){{M^ + Q^)/{M^ + Q^Y))-g{i,M^). 

Let us note that strictly speaking, we must use the generalized parton distributions (GPD) 



in Eq. 5.1 However the difference between gluon GPD and gluon pdf is not large in this case 
because fractions carried by gluons in GPD differ by the factor two at moderate x and 
tend to one at extremely large energies as the consequence of increase of parton momenta 
with energy. (In fact most of the non-diagonal effect in this approach is included in the 
wave functions of the initial and final photons.) Indeed, the DVCS amplitude is expressed 
in terms of gluon GPD where xi = M^/ s, X2 = (M^ + Q^) / s. Our analysis above shows that 
^ at moderate energies. So xi ~ x ~ 0:2/2 . Analysis of the QCD evolution equation 
for GPD [21] shows that xiG{Q^,xi,X2) = g{Q^, xi,X2) ^ g{Q^, {xi + X2)/2) = g{Q^, 1.5x). 
At small X g{Q^,3/2x)/g{Q^,x) ^ 1/(3/2)^/^ ^ 0.9. In the regime when xi ^ X2- 

So GPD coincides with the gluon density. As a result we may neglect the difference between 
GPD and distribution functions in the considered kinematics. The numerical analysis of 



Eq. 5.1 shows that the ratio R very slowly decreases with the increase of energy due to a 
slow increase of a ratio M^/Q^ discussed in the previous section, and R ~ 1.6 for HERA 
energies. 



The result Eq. 5J^ is however not complete since we neglected the contribution of the 
AJ configurations. Indeed, we already saw in the previos section (Table 3) that it gives 
significant contribution in a total cross-section. In the framework of the AJM model the ratio 
of amplitudes of the DVCS to DIS can be calculated within the leading twist approximation 
as [291: 

i?AjM = ^^log(l + g). (5.2) 
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Here the parameter ml = 0.3 — 0.5 GeV^ is the cut off parameter ttiq < m^, rrip is the p 
meson mass. 

Combining the pQCD and AJ model contributions we have 

_ -RpQCDCr + -Rajm<7ajm 2^ 
(Jt + ctajm 



Here the pQCD contribution into the total cross section is given by Eq. |2.7 and the 



contribution of AJ to the total cross section is given by AJM - Eq. |B1[ The results of 
numerical calculation with parameter mg ~ 0.5 GeV^ are given in Table 4. We depict these 
results as a function of x for several values of in Fig|4| The ratio R is close to 2 at HERA 
energies and increases with (from 5 to 100 GeV^ by ~ 40%). This result is in a good 
agreement with the analysis of the HI and ZEUS data in Ref. [31] (see in particular Table 
4 in Ref. PI]). Our main prediction is that the ratio R should decreases with the rise of 
energy. It tends to one at asymptotically large energies in agreement with the result for the 
BDR [23]. However the onset of this regime is very slow. This prediction can be checked 
experimentally in the study of DVCS processes at LHeC 

Our conclusion on the important role of AJM contribution in DVCS at HERA energies is 
in the qualitative agreement with the recent experimental data [32] that shows the important 
role of soft QCD in the diffractive processes in DIS at HERA. 

We want to draw attention that agreement between experimental results and theoretical 
prediction is rather good. This is due to the fact that the interaction of dipole effectively 
includes the NLO corrections since parton distributions were obtained by fitting the ex- 
perimental data. Consequently one may hope that NLO corrections to impact factors are 
relatively small. 

Let us stress that the current calculation is preliminary. More detailed calculation should 
account for the contribution of c-quark, and study in detail the dependence of R on the 
AJM parameters). 



VI. ON THE INCREASE WITH ENERGY OF PARTON MOMENTA IN THE 
CURRENT FRAGMENTATION REGION IN THE BLACK DISC LIMIT. 

In the previous sections we found the increase with the energy of the parton momenta 
in the current fragmentation region in the pQCD regime.. Here we shall consider the con- 
figurations in the projectile wave function that become dominant just near the transition 
to the black disk limit. We will explain that due to the dominance of these configurations, 
the increase with the energy of the parton transverse momenta in the current fragmentation 
region is further accelerated at the energies where the DGLAP and BFKL approximations 
to QCD violate conservation of probability/unitarity 

The aim of the present chapter is to analyze the properties of this new QCD regime (we 
also call it the black disk regime; let us stress that this regime is different from the possible 
new pQCD regime that was discussed in the previous sections). We shall study the increase 
of the characteristic transverse momenta of the configurations in the wave function of the 
projectile responsible for the blackening, and derive a expression for the total cross-section. 
This expression gives a total cross section as a function of a transverse momentum increase 
rate. 

We will use the probability conservation for the amplitude describing the dipole scattering 
at given impact parameter b. In the kinematics of a; <^ 1 Feynman amplitude describing scat- 
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tering of dipole off a target A{x, t, Q"^) is calculable in terms of partial amplitude x, Q'^) 
as 

A{x,t = -qlQ^) = 2sJ d^bexp{tqtb) f{b\x,Q^), (6.1) 

and 

fib^x,Q') = [ p^eM^m^^^^^. (6.2) 

J OTT"^ S 

In a quantum theory the elastic cross section for the scattering of dipole off target T is 
given by the well known formulae : 

aei = J d%\fib,x,Q')\\ (6.3) 

The inelastic cross section is given by: 



2, 



a,n = J d'b{l- l + tfib\x,Q') ), (6.4) 

and total cross sections is equal to: 

atot = 2 j d%lui f{b\x,Q^). (6.5) 

It follows from these equations that 

Im/(62,x,Q2) < 1. (6.6) 

The limit 

Im/(62,x,Q2) = i, (6.7) 

is reached in the small x regime when a small transverse size of the quark-gluon dipole is 
absorbed by the target with the 100% probability at given b. 

In the numerical evaluations it is convenient to use profile function Im/ = F and neglect 
Re/ as compared to /m/, which is justified at sufficiently high energies. 

These equations allow us to evaluate some important properties of the onset of the new 
QCD regime of strong interaction with small coupling constant and the kinematical region 
of applicability of perturbative QCD to hard processes. 



According to Eq,6^ the dependence of F on impact parameter b is given by the Fourie 
transform of the two gluon form factor of a nucleon. It was measured in the hard 
exclusive vector meson production at HERA and at FNAL. The two gluon formfactor 
can be parameterized as: F{q^) = 1/(1 + g^//i^)^, with /i decreasing with the decrease 
ofx (/i2(10-2) ~ 1 GeV^,fi\lO-^) ~ 0.7 GeV^), cf. discussion in Appendix C. Within 
the framework of pQCD at large kt the dependence of on a; due to Gribov diffusion 
in the impact parameter space should be a small effect. Variation of fi comes mainly 
from the diffusion at ~ Qq which is moderate, and it is further suppressed by the 
DGLAP evolution [31]. Thus 

r{b,x) = {xo/x)^^^^'^fibKi{fib). (6.8) 

Here Ki is the MacDonald function and kt is the parton transverse momentum within 
the dipole. Parameter Xo is function of kf calculable in pQCD. The dependence on en- 
ergy has been evaluated within pQCD or using the experimental data. Experimentally 
X{10GeV^) ^ 0.2 and increases with k^. 
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Knowledge of x, b, kt dependencies of T and Eq j6.7| allows us to evaluate the character- 
istic kt in the region, where the onset of the black disc regime occurs. This evaluation 
predicts significantly more rapid increase with energy of kt in the vicinity of black disc 
limit than the calculations performed in the preceding sections. This is the signature 
of the blackening of the interaction, i.e. of the transition to the regime described by 



Eq j6.7[ Since the integration over kt is logarithmic, it is determined by the upper limit 
of the integration over kf in the probability defined by the wave function of the dipole. 

The dependence of kt on collision energy is determined as the rate of increase of the 
black disk momenta squared with energy: 

kl ~ kl,{s/sor, (6.9) 

where Sq ~ 10'* GeV^ (corresponds to the HERA energy). The dependence of kfi, on s 
is depicted in Tables 5,6. Note that k does not depend on whether we take F = 1/2 
or r = 1 as a blackness condition, and is the same for fermionic and gluonic dipoles. 
Our results show that this coefficient weakly depends on energy, with k, = 0.3 ^ 0.4 at 
the LHC energies and beyond. 

The slope of t dependence of the cross sections of hard exclusive processes at t=0 - 
5(0) 

/ciW[^(6^x,Q^)] 

= IdWib,^,Q^)] ' ^'-''^ 
should increase with energy in the vicinity of BDR and at larger energies. This prop- 
erty follows from the increase with energy of the essential impact parameters. Using 



Eq,6^ we obtain equation for maximal b = bM where T is equal to 1: 

l/log(/i6Af^l(MM)) = Klog(xo/x). (6.11) 

At sufficiently large energies where k, log(a;o /x) > 1 the above equation has the solution: 

fibM = K\og{xo/x). (6.12) 

Evidently ^a/ (x — * 0) would become practically independent of xq. 

Note that within the LT approximation /j, and therefore the slope B are slowly in- 
creasing with energy in some interval of x ~ 0.01 as the consequence of spontaneously 
broken chiral symmetry. The long range pion cloud carries small fraction ^ niT^/mp 
of proton momentum and therefore scattering off this cloud is suppressed at not too 
small X because amplitude is rapidly decreasing with decrease of collision energy [I5] . 

A slightly more sophisticated reasoning which uses the completeness over the hadronic 



states shows that Eqs,6.11 6.12 are valid even if we account for the gluons and quark- 



antiquarks pairs in the wave function of the dipole. 



Knowledge of 6a/ and Eq j6.5| allows us to calculate the total cross section of the scat- 
tering of the spatially small dipole. We have F = 1 for 6 < 6a/, F(6 > 6a/) = 
6ii'i(/i(x)6)/6A/-ft'i(/u(x)6A/). Then we have 

a = 2nblj + 2 J d%Q{b - 6a/) F(6, x, Q^) 

= 277(6^^ + 2(6A///i)ir2(6M/i)/i^l(6M/i)). 

(6.13) 
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Here to calculate the second integral we use an identity between special functions: 



dx 



The second term in the second line of eq. 6.13 is the contribution of peripheral colli- 
sions. Thus cross sections of hard processes integrated over impact parameters would 
increase with energy forever. 

Note that using eq. 6.13 we can estimate, what part of the total cross section comes 
from the black disk regime, and what part from pQCD. Indeed, estimating Bm as a 
scale at which the black transverse momentum is of order ~ 1 GeV, we have 6m ~ 4.5 
GeV~l for typical x ~ 10-^ ~ 10 GeV^. We use a typical value of ~ 0.7 GeV 
for these energies. Then the contributions of black disk to total cross section is of 
the order 66%. Let us stress that here we assume Bm <^ l/(2m7r) which is a typical 
scale for the phenomenon of spontaneously broken chiral symmetry in the soft QCD 
physics. 



VII. THE COHERENCE LENGTH AND THE SHAPE OF THE FAST NUCLEON 
AND NUCLEI. 

A. The coherence length. 

In the previous sections we determined the energy dependence of the effective transverse 
scale at high energies. This result allows us to evaluate the coherence length. The coherence 
length Ic corresponds to the life-time of the dipole fluctuation at a given energy in the rest 
frame of the target. The original suggestion of the existence of the coherence length in the 
deep inelastic scattering was first made by loffe, Gribov and Pomeranchuk [211 ES]- It was 
found already in the sixties within the parton model approximation by loffe [36] that the 
coherence length at moderate a: is /c ~ l/2m]^x i.e. it linearly increases with energy. In the 
pQCD as a result of QCD evolution coherence length increases with energy more slowly: 

/, = (l/2m^x)(so/s)\ (7.1) 

Such a behavior of Ic with energy has been found previously in the numerical calculations 
of structure functions in the target rest frame [STJ [381 ES] • 

In a similar way, once we have the onset of the black limit regime, the coherence length 
corresponding to the dominant configurations in the process decreases even more slowly: 

h=(l/2m^x){s^/sY. (7.2) 

Here k is the rate of the increase of the black disk transverse momenta, given by eq. |6.9[ 



B. The shapes of the fast nucleons and nucleus. 

The results from the previous subsection have the important consequences for the trans- 
verse structure of the hadrons and nuclei. 

Let us consider the longitudinal distribution of the partons in a fast hadron. In the parton 
model the longitudinal spread of the gluonic cloud is ~ l//x for the wee partons (where 
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fi ~ 0.3 ^ 0.4GeV/c is the mass scale of soft QCD processes. ) and it is much larger than 
for harder partons, with ~ 1/xPh for partons carrying a finite x fraction of the hadron 
momentum pjOj . The picture is changed qualitatively in the limit of very high energies when 
interactions reach BD regime for kt ^ /i. In this case the smallest possible characteristic 
momenta in the frame where hadron is fast are of the order kt{BDR) which is a function 
of both initial energy and transverse coordinate, b of the hadron. Correspondingly, the 
longitudinal size is ~ l/kt{BDR) <C Note here that we are discussing longitudinal 

distribution for typical partons. There is always a tail to the momenta much smaller than 
typical all the way down to kt ^ fi which corresponds to the partons with much larger 
longitudinal size (a pancake of soft gluons corresponding to the Gribov's picture). However at 
large energies at the proximity of the unitarity limit the contribution of the gluons with kf < 
kth is strongly suppressed ^40j. In the BDR this tail is suppressed by a factor k^ /kt{BDR)^ 
in addition to the phase factor [2]. In the color glass condensate model the suppression is 
exponential [IT] . 

Since the gluon parton density decreases with the increase of b the longitudinal size of 
the hadron is larger for large b, so a hadron has a shape of biconcave lens, see Figs. 5][6 
In the numerical calculation we took 

|/,| = l/h{BDR), (7.3) 

neglecting overall factors of the order of one (typically in the Fourier transform one finds 
{z) ~ ^)- We calculated kt{BDR) for fixed external virtuality ~ iOGeV'^. Our results 

are not sensitive to the value of Qq, as the value of only enters in the combination 
x' = (Q^ + M'^)/s, and the k^ we found were comparable or larger than Q'^/4:. Indeed, 
the direct calculation shows that for small b the change of 1/kt if we go between external 
virtualities of 60 and 5 GeV^ is less than 5%. Such weak dependence continues almost to the 
boundary of the picture [5] where fc* ~ 1 GeV. Near the boundary the uncertainty increase to 
~ 25%, meaning that for large b (beyond those depicted in Fig. |5]) the nucleon once again 
becomes a pancake and there is a smooth transition between two pictures ( biconcave lens 
and pancake). We want to emphasize here that the discussed above weak dependence of 
kt{BDR) on the resolution scale indicates that the shape of the wave function for small x is 
almost insensitive to the scale of the probe. 

We depict the typical transverse quark structure of the fast nucleon in Fig. |5| We see that 
it is drastically different from the naive picture of a fast moving nucleon as a flat narrow 
disk with small constant thickness. (Similar plot for the gluon distribution is even more 
narrow). Note also that for the discussed small x range kt > IGeV/c for b < Ifm. Since 
the spontaneous chiral symmetry breaking corresponds to quark virtuality < IGeV"^, 



probably ~ 0.7Ge\^^ [12], corresponding to kt ~ y f/W^ ~ 0.7GeV^/c the chiral symmetry 
should be restored for a large range of b in the proton wave function for small x. 

Let us consider the case of the DIS on the nuclei. First, we consider the case of external 
virtualities of the order of several GeV. In this case the shadowing effects to the large extent 
cance the factor A^^^ in the gluon density of a nucleus for the central impact parameters, b 
[43] . and the gluon density in the nuclei is comparable to that in a single nucleon for 6 ~ 0. 
Consequently over the large range of the impact parameters the nucleus longitudinal size is 
approximately the same as in the nucleon at 6 ~ 0. 

However for very small x we find large kt{BDR) corresponding to Ak^{BDR) > 40 GeV^. 
This is a self consistent value as indeed for such the leading twist shadowing is small. 
Accordingly we calculated the shape of the nucleus for the external virtuality > 40 GeV^. 
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We should emphasize here that taking a smaller virtuality would not significantly change 
our result for kt{BDR) (at the same time LT nuclear shadowing reduces a low momentum 
tail of the kf distribution as compared to the nucleon case). In the discussed limit of the 
small leading twist shadowing, the corresponding gluon density unintegrated over b is given 
by a product of a nucleon gluon density and the nuclear profile function: 

T(6)= jdzp{h,z), (7.4) 

where the nuclear three-dimensional density is normalized to A. We use standard Fermi step 
parametrization |44] 

pir) = CiA) ^ Ra = I.IA^/Hm, a = 0.56fm. (7.5) 

^ ^ n + exp((r-i?A)/a) ^ ' ^ ^ 

Here r = + 6^, and A is the atomic number. C(A) is a normalization factor, that can 
be calculated numerically from the condition / d^rp{r) = A. At the zero impact parameter 
T(6) ^ 0.5^1/3 ig^^gg A. 

The dependence of the thickness of a fast nucleus as a function of the transverse size is 
depicted in Fig. |6] for a typical high energy s = 10'' GeV^, = 40 GeV^. We see that 
the nuclei also has a form of a biconcave lens instead of a fiat disk. The dependence on the 
external virtuality for the nuclei is qualitatively very similar to the case of the nucleon. For 
small b the dependence is very weak (of order 5%) and increases only close to the boundary 
of the biconcave lens region where it is of order 20% ( and A;^ ~ 1 GeV). For larger b we 
smoothly return to the pancake picture. 

Note that this picture is very counterintuitive: the thickness of a nucleus is smaller than 
of a nucleon in spite of ~ A^^^ nucleons at the same impact parameter. The resolution of the 
paradox in the BD regime is quite simple: the soft fields of individual nucleons destructively 
interfere cancelling each other. Besides for a given impact parameter b, the longitudinal size 
of a heavy nucleus l/k['^\BDR) < l/k^^\BDR) since the gluon distribution function in the 
nuclei Ga{,x^ b) > Gi\i{x, b). So a naive classical picture of a system build of the constituents 
being larger than each of the constituents is grossly violated. The higher density of partons 
leads to the restoration of the chiral symmetry in a broad b range and much larger x range 
than in the nucleon case. 

VIII. EXPERIMENTAL CONSEQUENCES. 

The current calculations of the cross sections of the hard processes at the LHC are based 
on the use of the DGLAP parton distributions and the application of the factorization 
theorem. Our results imply that the further analysis is needed to define the kinematic 
regions where one can use DGLAP distributions. We showed in the paper that for DIS at 
high energies in the kinematical region of sufficiently small x it is necessary to use the kt 
factorization and the dipole model, instead of the direct use of DGLAP. 

A similar analysis must be made for the pp collisions at LHC. It has been understood long 
ago that the probability of pp collisions at central impact parameter is close to 100% (total F 
is close to 1) even for soft QCD, i.e. at lesser energies than those necessary to achieve BDR for 
the hard interactions. The compatibility of probability conservation with the rapid increase 
of hard interactions with energy , predicted by QCD, requires the decrease of importance of 
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soft QCD contribution with energy |15]. As a result the hadronic state emerged in pp, pA, 
AA colhsions at sufficiently large energies consists of two phases. Central collisions would be 
dominated by the strong interaction with small coupling constant - the phase with unbroken 
chiral and conformal symmetries. On the contrary, the peripheral collisions are dominated 
by the more familiar phase with broken chiral and conformal symmetries. At these energies 
the QCD phase at central collisions - with the unbroken chiral and conformal symmetries 
-will be different from that for the peripheral collisions. This new phenomenona may appear 
especially important for the central heavy ion collisions at LHC and at RHIC. Quantitative 
analysis of this problem will be presented elsewhere. 

The hard processes initiated by the real photon can be directly observed in the ultraphe- 
ripheral collisions [16]. The processes where a real photon scatters on a target, and creates 
two jets with an invariant mass M^, can be analyzed in the dipole model by formally putting 

= 0, while is an invariant mass of the jets. In this case with a good accuracy the 
spectral density discussed above will give the spectrum of jets in the fragmentation region. 
Our results show that the jet distribution over the transverse momenta will be broad with 
the maximum moving towards larger transverse momenta with increase of the energy and 
centrality of the 7^4 collision. 

We have seen that our results can also describe DCVS processes. The ratio R of DCVS 
7* — i> 7* and forward amplitudes at t = is of order 2 at HERA energies at small external 
virtualities, and rapidly growing with Q^. This ratio slowly decreases with the decrease of 

X. 

Finally, our results can be checked directly, if and when the LHeC facility will be built 
at CERN. One of us, B.Blok, thanks S.Brodsky for the useful discussions of the results 
obtained in the paper. This work was supported in part by the US DOE Contract Number 
DE- FG02-93ER40771 and BSF. 

APPENDIX A: THE BASICS OF THE DIPOLE APPROXIMATION AP- 
PROACH. 

The QCD factorization theorem allows to express the total cross section of the DIS off a 
hadron through the convolution of the virtual photon wave function. LO approximation is 
equivalent to the dipole approximation and the cross section of the dipole scattering off a 
hadron [TU] . 

In this section we shall remind the reader the basic idea of the derivation of this result 
[ini ED]- Consider for the simplicity the case of the longitudinal photon. The total cross 
section of the photon scattering off the target T has the following factorized form: 



Here d^ji/k"^ is the propagator of the exchanged gluon in the light-cone gauge, TV is the sum 
of the box diagrams describing the scattering '~^*g qq, and T is the amplitude of the gluon 
scattering of the target. Restricting ourselves to the contribution of longitudinal gluon [17] 
polarization and using the Sudakoff decomposition, the latter equation can be rewritten as 




(Al) 



/ 



dadj3(Pkt 



{ImNf^p^pmimT^qf'qn. 



(A2) 



a = 



8(M)^(27r)4(A;2)2 
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Using the Ward identity we can rewrite: 



Since we integrate over d kt, the r.h.s. of Eq. A3 can be substituted by J2 1x^=1,2 ^^^jltvt^t 
The latter term can be directly connected with the scattering cross section — >■ qq: 

/it ft 

Here a, j3 are the Sudakoff parameters of the gluon momenta k: 

k = -aq' + pp + kt, (A5) 

q ia the photon and p is the target momenta. The total cross section then can be rewritten 
as (see ref ^20j for details) 

^r9jPGT{P,Q'). (A6) 

Thus the integrand in Eq. A2 for the total cross section is positively defined since it is a 
product of the cross section of the creation of qq pair and the structure function of a target. 
Here /5 = (Q^ + M^)/s. 

It is convenient to rewrite Eq. A6 by expressing the cross section cr{'y*g — qq) through 
the light-cone wave functions of the virtual photon: 



2(27r)3 J ' 2(27r)3 

X tp^rt, z){2tp{rt, z) - i){rt + kt) - tpin - h)) 

4:ImTli^^^^k'''k''^ 
s 

(A7) 

Here is the transverse momentum of the constituent within the dipole and kt is the 
transverse momentum of the exchanged gluon. The tensor T"^^ is the sum of the diagrams 
describing imaginary part of the amplitude for the gluon scattering off the target T. The 
integral of T"^^ over d'^kt is proportional to the gluon structure function of the target T. The 
function ip is the wave functions of the virtual photon . Here rt is transverse momentum 

Note that the validity of the dipole approximations is becoming better for the larger 
transverse momenta of the dipole constituents. 

Since rt ^ kt we can expand the wave functions in the Taylor series: 

2ij{rt, z) - ^{rt + kt) - ^{n - kt)) = -^{rt)Aij{rt)kl (A8) 



The latter equation is just the equation 2.1 



dz dhti'{rt)M{rt)g{{M^ + Q')/s, M^). (A9) 
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Here rf/{z{l — z)) = and z is the fraction of photon momentum carried by one of the 
constituents of the dipole. 

The resulting cross section can be written either in the momentum representation, giving 
the representation of the total cross section as the dispersion integral over the invariant 
masses of the dipoles, or in the coordinate representation, giving the cross section as the 
integral over the transverse dipole scales. 

Similar derivation can be made for the scattering of the spatially small transversely 
polarized photon. Principal difference is that the contribution of large size qq configuration 
is not suppressed in the case of transversely polarized photon -leading to the aligned jet 
model. 

Let us show now that within the leading logarithmic approximation one can rewrite Eq. 



A9 in the explicitly positive form. Indeed, let us perform the integration over d rt by parts 



using A = didi, i = 1, 2. The resulting expression contains four terms: 

a ^ I dzij d\di^pdiijgi{Q^ + M^)/sArt) + J d\d^^lj^Pirt)dig{iM^ + Q'')/sArt) 
+ boundary terms from — > + boundary terms, from rt ^ oo. 

(AlO) 

Here for the convenience we use definition: g = xG. Let us compare the relative contribu- 
tions of the second and first terms. The ratio of the corresponding integrands is determined 
by the ratio of derivatives of the logarithms of the wave function and of the structure func- 
tion. The differentiation {d/d{r^) removes one log(r^) from the pQCD series over the powers 
of as log(rj ) describing the pQCD evolution of parton distributions. As a result the logarith- 
mic derivative of parton distribution is suppressed by as- Hence in the leading logarithmic 
approximation the second term is parametrically smaller than the first one and can be ne- 
glected . Now consider the boundary term. The first of the boundary terms corresponds 
to the contribution of the small momenta (of — > 0), i.e. to the aligned jet model contri- 
bution. This contribution is known to be small for a^,. (It is zero within the parton model 
approximation). The second boundary term corresponds to the contribution of the very 
small dipoles, and clearly is zero since the tpdiipG ^ as ^ oo. We conclude, that in the 
leading logarithmic approximation only the first term should be retained, and we come to 



Eq. 2.3 in section 2 in the main text. 



APPENDIX B: THE ALIGNED JET MODEL (AJM) 

The purpose of this appendix is to remind a reader the basic formulae used to calculate 
the contribution of the AJM model to a total DIS cross sections and forward amplitude of 
DVCS process in the section 4. The contribution of the aligned jet model to the total cross 
section initiated by the transversely polarized virtual photon has been evaluated in |30j : 

= f "^^'p^2^-^-(^Ao)°-°^(^(M,^ - M^) + ^J{M^ - Ml)). (Bl) 

Here u is the upper cut off, Mq is a characteristic hadronic scale of the order IGeV^ . The 
first term in the formula corresponds to the contribution of the vector meson dominance 
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model and it is HT effect, while the second term is the aligned jet model contribution. Since 
our interest in this Appendix is focussed on moderately large we use R[M'^^ = p{e~^e~ 
hadrons)/p(e"''e~ — > fj,^fj,~)) = 2. The contribution of charm is not large but requires special 
treatment We assume that the cross section for the interaction of hadronic size qq pair with 
a nucleon target depends on collision energy as the soft Pomeron exchange, that the AJM 
and meson - nucleon cross sections are equal. The cross section of vector meson-proton 
interaction -avp is taken to be ~ 25 mb at the energies of the order sq ~ 400 GeV^. 

The contribution of the AJM model to the longitudinal model cross sections is suppressed 
by a factor 4 * u/Q^, i.e. 

4m 

0"LAJM ~ ^O^TAJM (B2) 



APPENDIX C: THE CHARACTERISTIC TRANSVERSE MOMENTA IN THE 
BLACK DISK LIMIT. 

Let us now evaluate the scale of the black disc limit where the partial wave for the forward 
scattering amplitude reaches 1. This scale gives us the biconcave shape of the nucleon, as it 
is explained in the section 5. 

The black disk scale was discussed already in a number of papers (see e.g. Ref. [16] for the 
recent discussion and references ) . Here we remind for a reader the rules of the calculation of 
the black scale in the dipole model, discuss it's dependence on energy. In order to calculate 
the corresponding we need to calculate the partial wave for a dipole -target scattering at 
zero impact parameter and then impose a black disc limit condition: the partial wave should 
be 1. (The scale of black disc limit depends on the impact parameter, so in the analysis 
below we shall consider the case of the zero impact parameter, for larger impact parameters 
the black disc limit arises at higher energies (see section 6). 

Let us define in the standard way the profile function for the dipole scattering: 

r(s,6) = / d\exp{tqb)A{s,t) (CI) 

Here t = —qf, and A(s,t) is the amplitude of the dipole scattering on the nucleon with the 
momentum transfer q. The elastic, inelastic and total crosssections are connected with this 
profile function as 

atot(s) = 2j d^bReT{s,b) 
aei(s) = J d%\r{s,b)\' 
(Tin(s) = J d%{2Rer{s,b) -\r{s,b)\^). 

(C2) 

The amplitude of dipole nucleon scattering can be expressed through the total cross 
section and gluon density 

A{s,t) =isa{s)Fg{t). (C3) 
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Here a{s) is given by Eq. 1.2 with = 4/3 for the fermion dipole and = 9/4 for the 
gluonic dipole. The two-gluon form factor of the nucleon can be parametrized as to satisfy 
mainstream ideas on nucleon form factors: 



For the numerical dependence of the two-gluon nucleon form factor parameter we 
shall use the parametrization 



ml{x,M') = ml{x,Ql){l + 1.51og(gVg^))°-°°9^°^(i/-), (C5) 



where 



ml{x,Ql) = 8/C,C = max(0.28fm^0.031/m2 + 0.0194 log(0.1/x)fm2), (C6) 

and l/m2=0.04 GeV-^. Qualitatively, ~ 1 GeV^. 

This parametrization is based on the fit of cross sections of hard diffractive processes in 
the kinematical domain: ~ 10^^ — 10^^ and Q^up to 10^ GeV^ made in ref. [T3] and was 
used in ref. |16] for the analysis of LHC collisions. 

Since absorption can not exceed 100% : 

r(s,6)<l (C7) 



Let us make the Fourier transform of the form factor C4 into impact parameter space : 



F{b) = {l/2nn J cPqeMmH/H + q',f = (l/2vr)^^^^^^^^ (C8) 
Then the condition of probability conservation (of 100% absorption) can be rewritten as 

r(s,i) = '^^M,n,b)<l (C9) 

Note that described above approach to the probability conservation differs from the satu- 
ration models which explore formulae of LO BFKL approximation derived for hard processes 
with one scale and accounted for the conservation of probability in terms of elastic eikonal 
approximation. ( For the recent review see ^81 H9] ) .Our approach is based on the DGLAP 
formulae derived for two scale (hard-soft scales) processes like total cross section of DIS 
and on the condition of complete absorption whose validity does not require applicability 
of eikonal approximation. Note also that in QCD inelastic diffraction in the scattering of 
spatially small dipole is significantly larger than elastic one [15]. For many phenomena (but 
not for all) both approaches lead to qualitatively similar but quantitatively different pre- 
dictions . Another advantage of the approach used in the paper is in the exploring impact 
parameter gluon distribution which follows from the two gluon form factor measured in the 
hard diffractive processes [15]. 

The typical dependence of the BDR onset scale fc^ on energy is presented in Table 5 
for the gluonic dipole scattering off a proton target, and in table 6 for fermionic dipole, 
where we used a CTEQ5 structure functions, the use of CTEQ6 leads to similar results. 
We consider two possible criterions for reaching the black disk limit: one when the partial 
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wave r at the central impact parameter reaches 1, another when it reaches 1/2 [15]. Indeed, 
when the partial wave reaches 1/2 the probability of inelastic interactions reaches 3/4, i.e. 
interactions become strong and pQCD can not be used any more \15\ . 

Let us note that the interaction of quark dipoles (see table 6) is far from the black disk 
regime for HERA energies, although nonperturbative effects (corresponding to F = 1/2) 
seem to start to appear for kt < 1 GeV. The results in tables 5,6 are in good agreement with 
the previous determination of these scales fiSl [T5l H6] . 

For our purposes it will be important to determine the rate of increase of the fc^. This 
rate does not depend on the type of the dipole (fermionic or gluonic) or on the partial wave 
condition. 

For the gluonic dipole in DIS for s up to s ~ lO'' GeV^ the BDR scales increase relatively 
slow as for CTEQ6 and s"'^ for CTEQ5 parametrisation. However at the energies 
above s ~ 10^ GeV^ the increase speeds up with s^,n = 0.35 starting from s ~ 

10^ GeV^ energies , n=0.32 for the exponent for CTEQ6 parametrisation. The CTEQ5 
parametrisation leads to even more rapid increase with the energy with n ~ 0.4. 

APPENDIX D: THE PARAMETRISATION OF THE GLUON DISTRIBUTION 
FUNCTION 

Throughout this paper we use the following approximate formulae for CTEQ5L and 
CTEQ6L distribution functions. 

The CTEQ5L distributions has been shown to be in good agreement with HERA data 
|16j . in the range of x ~ 10~^, 10~^. To extrapolate to very small x, we shall use the 
approximate formulae in the form: 

^~a(M2)/x^(*^'). (Dl) 

Here the functions 

a{M^) = 2.00123 - 1.69772 ■ 10/M^ + 3.07651/yMVia - 0.228087 ■ log (M^/IO), (D2) 

c(M2) = 0.045 \og{M^) + 0.17, (D3) 

where is in GeV^. This formula is also the good fit to the observed behavior of the 
structure functions measured at HERA for 150GeV^ > > 3 GeV^ , made by ZEUS and 
HI collaborations [50] . 

For CTEQ6 we have similar approximation, but with 

a{M^) = 3.96 - 0.032 ■ lO/M^ - 1.82/^MVlO. - 0.527 ■ log (M710), (D4) 

c{M^) = 0.035 \og{M^) + 0.2, (D5) 

The CTEQ6 distribution is qualitatively similar to CTEQ5 but grows slightly more slow, 
and is in better agreement with the J\l/ production data. 
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APPENDIX E: APPROXIMATE FORMULAE FOR THE MEDIAN TRANS- 
VERSE MOMENTA. 



For longitudinal photons it is easy to see that the results in a Table 1 with a very good 
accuracy can be described by an approximate formula 



2 



(x/0.01)°-°^+°-°2i°s('3VQo- ' 
a(g2) = 0.6GeV2 + 0.ig2. 

(El) 

Here = 10 GeV^. 

For the transverse photons the analysis of Table 2 gives an approximate formula 



(x/0.01)°-°^+°-°^^'°s('5^/'5o) ' 
a(Q2) = O.SGeV^ + 0.025Q^ 

(E2) 

The invariant masses increase slightly more slowly, with M'^{x ~ 0.01) ~ 0.95Q^ + 9GeV^, 
and M2 ~ i/(a./o.01)°-°^+°-°^5i°s(QVQg)^ Here Ql = 5 GeV^. If we use the CTEQ5 we 



obtain qualitatively the same results as in Eq. E2 but the exponent is slightly larger 

i2 /r)2^ 



Akf ~ /Qo)^ and the invariant masses increase also more rapidly, like 

(0.95Q2 + 9GeV2))/(a;/0.01)008+ooi5iog(QVQg). 
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g2 = 3 GeV^ 


5 GeV^ 


10 GeV^ 


20 GeV^ 


40 GeV^ 


60 GeV^ 


80 GeV^ 


100 GeV^ 


s = 10-^ GeV^ 


0.82 GeV^ 


1.1 GeV^ 


1.8 GeV^ 


3.14 GeV^ 


5.5 GeV^ 


7.75 GeV^ 


9.9 GeV^ 


11.9 GeV^ 


s = 10-4 GeV2 


0.90 GeV^ 


1.25 GeV2 


2.1 GeV2 


3.65 GeV^ 


6.5 GeV2 


9.3 GeV2 


11.9 GeV^ 


14.5 GeV2 


s = 10-5 GeV2 


0.95 GeV2 


1.35 GeV2 


2.3 GeV2 


4 GeV2 


7.2 GcV2 


10.2 GcV2 


13.15 GeV2 


16 GcV2 


s = 10-6 GeV2 


1.01 GeV2 


1.46 GeV2 


2.5 GeV2 


4.4 GeV2 


7.9 GeV2 


11.2 GeV2 


14.5 GeV2 


17.6 GeV2 


s = 10-7 GeV2 


1.08 GeV2 


1.61 GeV2 


2.7 GeV2 


4.8 GeV2 


8.65 GeV2 


12.3 GeV2 


15.9 GeV2 


19.4 GeV2 



TABLE I: The scale k^^ (50% of the total cross section) for longitudinal photons in DIS (CTEQ6) 





Q2 = 3 GeV2 


5 GeV2 


10 GeV2 


20 GeV2 


40 GeV2 


60 GeV2 


80 GeV^ 


100 GeV^ 


s = 10-^ GeV2 


0.91 GeV^ 


1.07 GeV2 


1.35 GeV2 


1.78 GeV2 


2.3 GeV2 


2.75 GeV2 


3.06 GeV2 


3.45 GeV^ 


s = 10-4 GeV2 


1.18 GeV^ 


1.42 GeV2 


1.82 GeV2 


2.42 GeV^ 


3.25 GeV2 


3.9 GeV2 


4.4 GeV^ 


5 GeV2 


s = 10-5 GcV^ 


1.37 GeV2 


1.67 GeV2 


2.18 GeV^ 


2.98 GeV2 


4.1 GeV2 


5. GeV2 


5.8 GeV2 


6.5 GeV2 


s = 10-6 GcV^ 


1.54 GeV^ 


1.91 GcV^ 


2.57 GeV^ 


3.57 GeV2 


5. GeV2 


6.3 GeV2 


7.3 GcV2 


8.3 GeV2 


s = 10-7 GeV2 


1.76GeV2 


2.21 GeV^ 


2.95 GeV^ 


4.3 GeV2 


6.2 GeV2 


7.9 GeV2 


9.35 GeV^ 


10.5 GeV^ 



TABLE IL The scale k^^ (50% of the total cross section) for transverse photons in DIS (CTEQ6) 





Q2 = 3 GeV2 


5 GeV2 


10 GeV2 


20 GeV2 


40 GeV2 


60 GeV^ 


80 GeV^ 


100 GeV^ 


X = 10-2 GeV2 


38 


45 


50 


56 


60 


61 


62.5 


63 


X = W-^ GcV^ 


61 


66 


71 


75 


78 


80 


80.5 


81 


X = 10-4 GcV2 


73 


81 


84 


84.5 


87 


88 


88.5 


89 


X = 10-5 GeV2 


81 


84 


87.5 


90 


92 


92.5 


93 


93.5 


X = 10-6 GeV^ 


87 


89.5 


92 


94 


95 


95.5 


96 


96 


X = 10-7 GeV2 


91 


93 


95 


96 


97 


97.3 


97.6 


98 



TABLE III: The pQCD contribution (%) in the total cross section, that includes AJM contribution 
(CTEQ6) 
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X = 10"4 GeV2 


1.74 


1.82 


1.88 


1.92 


1.94 


1.95 


1.97 


X = 10-5 GeV2 


1.66 


1.72 


1.76 


1.79 


1.8 


1.81 


1.82 


X = 10-6 GeV^ 


1.6 


1.63 


1.67 


1.69 


1.71 


1.72 


1.72 


X = IQ-^ GeV^ 


1.52 


1.56 


1.6 


1.62 


1.64 


1.66 


1.66 



TABLE IV: The DVCS ratio R (CTEQ6) 





r = 1 


r = 1/2 


s = 10'^ GcV2 


1.6 GeV2 


3.1 GeV2 


s = 10^ GeV2 


3 GcV2 


5.6 GcV2 


s = 10^ GeV2 


6 GcV2 


10 GcV2 


s = 10^ GeV2 


11 GeV^ 


21 GeV2 



TABLE V: The scale for the onset of the black disk regime for the gluon dipole. 





r = 1 


r = 1/2 


s = 10^ GeV^ 


-GeV2 


1 GeV^ 


s = 10^ GeV^ 


-GeV2 


2 GeV^ 


s = 10^ GeV^ 


2 GeV^ 


5 GeV^ 


s = 10^ GeV2 


5 GeV2 


10 GeV2 



TABLE VL The scale for the onset of the black disk regime for the fermionic dipole. 
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FIG. 3: The contribution R of pQCD to the total cross section, that is a sum of pQCD and AJM 
model contributions. The cut off of the AJM model is 0.35 GeV^, for = 5 GeV^ (lower curve), 
20 GeV^ (middle curve), 40 GeV^, the x axis corresponds to log;^o(l/2^) 
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FIG. 5: 3D image of the fast nucleoli at s = 10^ GeV^ and the resolution scale < 40GeV^, 
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FIG. 6: 3D image of the fast heavy nucleus (gold) at s = 10^ GeV^ and the resolution scale 
q2 ^ 40GeV^.Let us stress that the thickness of a nuclei in the Figure is much smaller than the 
thickness of the nucleon. 
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